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Iterated functions system (IFS) is defined by specifying a set of functions in a classical phase 
space, which act randomly on an initial point. In an analogous way, we define a quantum iterated 
functions system (QIFS), where functions act randomly with prescribed probabilities in the Hilbert 
space. In a more general setting a QIFS consists of completely positive maps acting in the space of 
density operators. This formalism is designed to describe certain problems of nonunitary quantum 
dynamics. We present exemplary classical IFSs, the invariant measure of which exhibits fractal 
structure, and study properties of the corresponding QIFSs and their invariant states. 
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I. INTRODUCTION 



An iterated function system (IFS) may be considered as a generalization of a classical dynamical system, which 
permits a certain degree of stochasticity. It is defined by a set of k functions fi : fl ^ il, i = 1, . . . , k, which represent 
discrete dynamical systems in the classical phase space fl. The functions fi act randomly with given place-dependent 
probabilities pi : Q ^ [0,1], i = 1, ... ,k, X^iLi Pi = 1 Q- They characterize the likelihood of choosing a particular 
map at each step of the time evolution of the system. 

There exist different ways of investigating such random systems. Having defined an IFS, one may ask, how is an 
initial point xo & ^ transformed by the random system. In a more general approach, one may pose a question that 
how does a probability measure fi on fl change under the action of the Markov operator P associated with the IFS. If 
the phase space is compact, the functions fi are strongly contracting, and the probabilities pi are Holder continuous 
and positive (i.e. pi > 0), then there exists a unique invariant measure fi* of P - see for instance 0, (2, 1^], and 
references therein. 

For a large class of IFSs, the invariant measure /i* has a fractal structure. Such IFSs may be used to generate 
fractal sets in the space SI. In particular, iterated function systems leading to well-known fractal sets, such as the 
Cantor set or the Sierpiiiski gasket, can be found in Refjl|[. These intriguing properties of IFSs allowed one to apply 
them for image compression, processing, and encoding 00,0- 

Iterated function systems can also be used to describe several physical problems, where deterministic dynamics 
is combined with the random choice of interaction. In particular, IFSs belong to a larger class of random systems 
studied in Ref . 0, Q| . Such a composition of deterministic and stochastic behavior is important in numerous fields of 
science, since very often an investigated dynamical system is subjected to an external noise. 

Nondeterministic dynamics may also be relevant from the point of view of quantum mechanics. Although unitary 
time evolution of a closed quantum system is purely deterministic, the problem changes if one tries to take into 
account processes of quantum measurement or a possible coupling with a classical system. In the approach of Event 
Enhanced Quantum Theory (EEQT) developed by Blanchard and Jadczyk Q, the quantum time evolution is piecewise 
deterministic and in certain cases may be put into the framework of iterated function systems While some 
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recent investigations in this area concentrate mostly on IFS's acting in the space of pure states we advocate a 
more general setup, in which IFS's act in the space of mixed quantum states. 

The main objective of this paper is to propose a general definition of quantum iterated function system (QIFSs). 
Formally, it suffices to consider the standard definition of IFS and to take for fl an iV-dimensional Hilbert space Hn- 
Instead of functions fi, i = 1, . . . , A:, representing classical maps, one should use linear functions Fi : Hn ^ Ti-N , which 
represent the corresponding quantum maps. Alternatively, one may consider the space Mn of density matrices of size 
N and construct an iterated function system out of k positive maps Gi : Mn ^ -Mn ■ The QIFSs defined in this way 
can be used to describe processes of quantum measurements, decoherence, and dissipation. Moreover, QIFSs offer 
an attractive field of research on the semiclassical limit of quantum random systems. In particular, it is interesting 
to explore quantum analogues of classical IFSs, which lead to fractal invariant measures, and to investigate, how do 
quantum effects smear fractal structures out. 

This paper is organized as follows. In the following section we recall the definition and basic properties of the classical 
IFSs, and discuss several examples. In Sec. Ill we propose the definition of QIFSs, investigate their properties, and 
relate them to the notion of quantum channels and complete positive maps used in the theory of quantum dynamical 
semigroups. The quantum-classical correspondence is a subject of Sec. IV, in which we compare dynamics of exemplary 
IFSs and the related QIFSs. Concluding remarks are presented in Sec. V. 

II. CLASSICAL ITERATED FUNCTION SYSTEMS 

Consider a compact metric space fl and k functions fi : ^ fl, where i — 1, . . . ,k. Let us specify k probability 
functions pi : Q ^ [0, 1] such that for each point x £ fl the condition X^iLi Pii^) = 1 is fulfilled. Then the functions fi 
may be regarded as classical maps, which act randomly with probabilities pi. The set Tci ■= {^i fiiPi '■ i = ■ ■ ■ , k} 
is called an iterated function system (IFS). 

Let M{fl) denotes the space of all probability measures on Q. The IFS !Fci generates the following Markov operator 
P acting on M (O) 

{Pfi)iB)^J2[ Prix)Mx) , (1) 

where B is a measurable subset of fl and a measure /i belongs to M{fl). This operator represents the corresponding 
Markov stochastic process defined on the code space consisting of infinite sequences built out of k letters which label 
each map fi. On the other hand, P describes the evolution of prohahility measures under the action of Tci- 

Consider an IFS defined on an interval in E and consisting of invertible maps {fi : i = 1, . . . , fc}. This IFS 
generates the associated Markov operator P on the space of densities ,I2||, which describes one step evolution of a 
classical density 7 

i 

where for x Q the sum goes over i = 1, . . . , k, such that x € fii^)- 

Let d{x, y) denotes the distance between two points x and y in the metric space fi. An IFS Tg\ is called hyperbolic, 
if it fulfills the following conditions for all i = 1, fc: 

(i) fi are Lipschitz functions with the Lipschitz constants Li < 1, i.e., they satisfy the contraction condition 
d{fi{x),fi{y)) < Lid{x,y) for all x,y Gil; 

(ii) the probabilities pi are Holder continuous, i.e., they fulfill the condition \pi (x) — pi (y)| < Kid{x,y)" for some 
a e (0, 1], K^e R+ for all x,yen; 

(iii) all probabilities are positive, i.e., Pi{x) > for any a; G il. 

The Markov operator P associated with a hyperbolic IFS has a unique invariant probability measure /z* satisfying 
the equation P/x, = /x*. This measure is attractive, i.e., P"/i converges weakly to /i* for every /i £ M{fl) as n ^ 00. 
In other words, u dP"^ tends to u d/i* for every continuous function u : M. Let us mention that the 

hyperbolicity conditions (i)-(iii) are not necessary to assure the existence of a unique invariant probability measure - 
some other, less restrictive, sufficient assumptions were analyzed in Refs. 14, 15, 16, 17). 

Observe that in the above case, in order to obtain the exact value of an integral /^^ u dfi^, it is sufficient to find the 
limit of the sequence u d{P^^) for an arbitrary initial measure /i. This method of computing integrals over the 
invariant measure /z* is purely deterministic Q. Sometimes it is possible to perform the integration over the invariant 
measure analytically, even though /x* displays fractal properties [l8] . Alternatively a random iterated algorithm may 
be employed by generating a random sequence Xj € by the IFS, j = 0, 1, . . . , which originates from an arbitrary 
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initial point xq. Due to the ergodic theorem for IFSs [3. Hfll \2(\. the mean value (!/«) X]j=o^ "(^j) converges with 
probability one in the limit n — > oo to the desired integral u dfi^ for a large class of u. 

If probabilities Pi are constant we say that an IPS is of the first kind. Such IFSs are often studied in the mathematical 
literature (see Ref. 0| and references therein). Moreover they have also some applications in physics. For example, 
they were used to construct multifractal energy spectra of certain quantum systems [23|, and to investigate second 
order phase transitions [22]. On the other hand, IFSs with place-dependent probabilities can be associated with some 
classical and quantum dynamical systems 0, US 0, EM EE E3| ■ In analogy with the position-dependent gauge 
transformations such IFSs may be called iterated function systems of the second kind • 

If 17 is a compact subset of R", while dE{x,y) represents the Euclidean distance, or is a compact manifold (e.g. 
sphere 5^ or torus T") equipped with the natural (Riemannian) distance dR, then an IFS will be called classical. 
For concreteness we provide below some examples of classical IFSs. The first example demonstrates that even simple 
linear maps fi may lead to a nontrivial structure of the invariant measure. 

Example 1. = [0, 1], fc = 2, pi = p2 = 1/2 and two afHne transformations are given by fi{x) — x/3 and f2{x) = 
x/3 -I- 2/3 for a; e fi. Since both functions are continuous contractions with Lipschitz constants Li = ^2.= 1/3 < 1, 
this IFS is hyperbolic. Thus, there exists a unique attracting invariant measure /j,*. It is easy to show |l| that /i* is 
concentrated uniformly on the Cantor set of the fractal dimension d = In 2/ In 3. 

The next example presents an IFS of the second kind. 

Example 2. As before, fl = [0, 1], k ~2, ,fi{x) = x/3, and f2{x) = x/3 + 2/3 for x fl. The probabilities are now 
place dependent, pi{x) = x and P2{x) — 1 ~ x. Although this IFS is not hyperbolic (condition (iii) is not fulfilled) , a 
unique invariant measure /i* still exists. It is also concentrated on the Cantor set, but now in a non-uniform way [l^ . 
The measure /i* displays in this case multifractal properties, since its generalized dimension depends on the Renyi 
parameter. 

Example 3. = [0, 1] x [0,1] C M.^ , k — 4, pi ~ p2 = ps = Pi ^ l/i. Four affine transformations are given by 




Also, this IFS is not hyperbolic, since the transformations /, are not globally contracting, the former two contract 
along X-axis, while the latter two along y axis only. An invariant measure /z* for this IFS is presented in Fig.QJl. The 
support of /X* is the Cartesian product of two Cantor sets. Thus, its fractal dimension is d= 2 In 2/ In 3. 

Example 4. Let = S*^. Take k = 2, pi = p2 = 1/2, and choose /i to be the rotation along z-axis by angle 
Xi (later referred to as Rz{xi))- In the standard spherical coordinates, /i(0, 0) = {9,(f> + xi)- The second function 
/2 is a rotation by angle X2 along an axis inclined by angle /? with respect to z-axis. Since both classical maps are 
isometries this IFS is by no means hyperbolic. The properties of the Markov operator depend on the angle (3, and 
the commensurability of the angles Xi- However, the Lebesgue measure on the sphere is always an invariant measure 
for this IFS. 

Example 5. 17 = [0, 1], fc = 2, pi ^ p2 = 1/2, Mx) = 2x for x < 1/2, and fi{x) = 2(1 - x) for a; > 1/2 {tent 
map); f2{x) = 2x for x < 1/2 and f2{x) = 2a; — 1 for a; > 1/2 {BernouUi map). Both classical maps are expanding 
(and chaotic), thus the IFS is not hyperbolic. The Lebesgue measure in [0, 1] is an invariant measure /i* for this IFS. 



III. QUANTUM ITERATED FUNCTION SYSTEMS 
A. Pure states QIFSs 

To describe a quantum dynamical system we consider a complex Hilbert space Ti.. When the corresponding classical 
phase space is compact, the Hilbert space Hn 'is finite dimensional and its dimension N is inversely proportional 
to the Planck constant h measured in the units of the volume of ft. Analyzing quantum systems, N is usually treated 
as a free parameter, and the semiclassical limit is studied by letting N ^ oo. 

A quantum state can be described by an element \ip) oi Hn normalized according to {ip\ip) = 1. Since for any 
phase a the element \^p') = e*"|V') describes the same physical state as we identify them, and so the space of all 
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pure states Vn has 2N — 2 real dimensions. From the topological point of view, it can be represented as the complex 
projective space CP^^^ equipped with the Fubini-Study (FS) metric given by 

i?F5 (1 0), I V') ) =arccos 1(01 V') I . (4) 

It varies from zero for \4>) = j-i/') to 7r/2 for any two orthogonal states. In the simplest case of a two-dimensional 
Hilbert space Hi the space of pure states Vi reduces to the Bloch sphere, CP^ ~ 5^, and the FS distance between 
any two quantum states equals to the natural (Riemannian) distance between the corresponding points on the sphere 
of radius 1/2. 

Definition 1. To define a (pure states) quantum iterated function system (QIFS) it is sufficient to use the general 
definition of IFS given in Sect. II, taking for fl the space Vn- We specify two sets of k linear invertible operators: 

• Vi : Hn Ti-N {i — ^, ■ ■ ■ ,k), which generates maps Fi : Vn Vn {i — ^, ■ ■ ■ , k) hy 

• Wi : Hn — > Ti-N (i = 1, ■ • • , fc), forming an operational resolution of identity, X^iLi W}Wi = 1, which generates 
probabilities pi : Vn ~> [0, 1] (i = 1, • ■ ■ , fc) by 

pdm--=\mimf (6) 

for any \(p) G Vn- 

Clearly, for any |0) g Vn the normalization condition = 1 is fulfilled. In this situation a QIFS may 

be defined as a set 

Tn^{Vn; F,:Vn^Vn; p^:PJv^[0,l]:^ = l,...,fc}. (7) 



Such a QIFS may be realized by choosing an initial state |0o) G Vn and generating randomly a sequence of pure 
states (|0j))jgN. The state l^o) is transformed into |0i) = Pi(|0o)) with probability pi(|(/)o)), later is mapped into 
\4>2) = Fj with probability pj (|0i)), and so on. If we choose Wi = ^/pil, then the probabilities are constant: 

Pi — Pi for i = 1, - - - , k. An arbitrary QIFS J-'n determines by formula ^ the operator P acting on probability 
measures on Vn- 

Such defined QIFS J^n cannot be hyperboHc, since the quantum map Fi are not contractions with respect to the 
Fubini-Study distance in cuIPn- 

Example 6. fl ^ Vn ^ CP^"^, fc = 2, pi = = 1/2, Fi{\tP)) = Ui{\ip)) and FiQip)) = UiQ^)), where the 
operators Ui {i = 1,2) are unitary. In this case both quantum maps are isometrics. Thus the natural Riemannian 
(Fubini-Study) measure in Vn is invariant, but as we shall see in the next section, its uniqueness depends on the 
choice of t/i and U2 - 



B. Mixed states QIFSs 

Mixed states are described by TV— dimensional density operators p, i.e., positive Hermitian operators acting in H^/ 
with trace normalized to unity, p — p\ p > and trp — 1. They may be represented (in a non unique way) as a 
convex combination of projectors. We shall denote the space of density operators by A^at. 

Definition 2. Now we can formulate the general definition of a QIFS as a set 

^N {Mn, G,:Mn^ Mn, p^ ■- Mn ^ [0, 1]; i - 1, k} , (8) 

where the maps Gi, i = 1, . . . , fc transform density operators into density operators, and for every density operator 
p G Mn the probabilities are normalized, i.e., J2i=iPiip) — ^■ 

The above definition of QIFS is more general than the previous one, since in particular d and pi may be defined 

by 

G.ip) = ^^ (9) 
tr (v^pV^j 
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and 

(p) = tr (w.pW^^ (10) 

for i — 1, . . . , k and p G Mn, where the linear maps Vi and Wi are as in Definition 1. Thus, each QIFS on Vn can 
be extended to a QIFS on A^at. Note that in this case pi (p) = triW^Wip). Hence, we can alternatively define the 
probabilities by Pi{p) = tr (Cip) {i ~ 1, . . . ,k, p £ Mn), where the linear operators Ci are Hermitian, positive, and 
fulfill the identity X^iLi — ^■ 

Now the dynamics takes place in the convex body of all density matrices Mn- 
The space of mixed states A4n has N'^ — 1 real dimensions in contrast to the 
(27V — 2)— dimensional space of pure states Vn- For = 2 its is just the 3-dimensional Bloch ball, i.e., the 
volume bounded by the Bloch sphere. 

The special class of QIFSs is a class of homogenous QIFSs introduced in more general setting by one of the authors 
[2^ . A QIFS is called homogenous if both pi and Gi -pi are affine maps for i = 1, . . . ,k. The mixed states QIFS being 
a generalization of a pure state QIFS and defined by formulas and lfT?Hl is homogenous if = for i = 1, . . . , fc. 
Interesting examples of such systems acting on the Bloch sphere where recently analyzed by Jadczyk and Oberg 
For a homogenous QIFS Pi and Gi may be interpreted in terms of a discrete measurement process as the probability 
that the measurement outcome is i, and the state of the system after the measurement if the result was actually i, 
respectively. 

A homogenous QIFS generates not only the Markov operator P acting in the space of probability measures on A^at, 
but also the linear, trace-preserving, and positive operator A: Mn ^ Mn defined by 

fc k 

Hp) = Y.P^{p)G,{p) = ^^P^^ (11) 

i=l 1=1 

iov p £ Mn ■ 

A mixed state p is A— invariant if and only if it is the barycenter of some P— invariant measure p, i.e., 

p = [ pdp{p) , (12) 

see Ref. 

Example 7. ft — Mn , k — 2, pi — p2 = 1/2, Gi{p) = UipUl and G2{p) — U2pUl. This is just Example 6 in 
other casting; the normalized identity matrix, — \/N is A— invariant irrespectively of the form of unitary operators 
Ui,i = 1,2. Note that p = p* may be represented as Eq. (|12ll . where the measure Ji, uniformly spread over Vn (the 
Fabini-Study measure), is P invariant. 

To define hyperbolic QIFSs one needs to specify a distance in the space of mixed quantum states. There exist 
several different metrics in Mn, which may be applicable (see e.g. Ref. |j2^|22i and references therein). The standard 
distances: the Hilbert-Schmidt distance 

Dhs{puP2) = Vtr[(pi - P2)2] , (13) 

the trace distance 

Ar(/0l,P2)=trV(pi-p2)2 -IIP1-P2 11^^ , (14) 

and the Bures distance 

DBures (pi , P2) - ^2 { 1 - tr [(p^ ' ) ^/^l } (15) 

the latter based on the idea of purification of mixed quantum states ^3r,'3^, are mutually bounded 113 • They generate 
the same natural topology in A^jv. Having endowed the space of mixed state with a metric, we may formulate 
immediate conclusion from the theorem on hyperbolic IFSs. We define a hyperbolic QIFS as in the previous section, 
and the following proposition holds. 

Proposition 1. If a QIFS 0) is homogenous and hyperbolic (that is, the quantum maps d are contractions with 
respect to one of the standard distances in Mn, Pi are Holder continuous and positive), then the associated Markov 
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operator P possesses a unique invariant measure ^. This invariant measure determines a unique A— invariant mixed 
state p & Mn given by Eq. 

Note that for a homogenous hyperboHc QIFS, the sequence A"(po) tends in the Hmit n ^ oo to a unique invariant 
state p irrespectively of the choice of an initial state po • 

Example 8. Q. = Aijq , k = 2, pi = p2 = 1/2, Gi(p) = (p + 2pi)/3 and G2{p) = {p + 2p2)/3, where we choose 
the both projectors pi — |1)(1| and pi — |2)(2| to be orthogonal. Since both homotheties Gi are contractions (with 
the Lipschitz constants 1/3) this QIFS is hyperboHc and a unique invariant measure ]! exists. In analogy with the 
IFS discussed in Example 1 we see that the support of Jl covers the Cantor set at the line joining both projectors pi 
and p2. However, this is nothing but a rather sophisticated representation of the maximally mixed two- level state 
P* — (Pi + P2)/2, which follows from the symmetry of the Cantor set and may be formally verified by performing the 
integration prescribed by Eq. ltT2|l . 



C. Completely positive maps and unitary QIFSs 

From the mathematical point of view it may be sufHcient to require that the map A is positive, that is, it transforms 
a positive operator into another positive operator. From the physical point of view it is desirable to require a stronger 
condition of complete positivity related to a possible coupling of the quantum system under consideration with an 
environment. A map A is completely positive (CP-map), if the extended map A (g) 1 is positive for any extension of 
the initial Hilbert space, TCn Ti-N '^'He, which describes coupling to the environment j33.l33|. 

It is well known that each trace preserving CP-map A (sometimes called quantum channel), can be represented 
(non uniquely) in the following Stinespring-Kraus form 

k k 

p' = Kk{p)^Y.'^jPVI with ^1/.V, = 1, (16) 

where Hnear operators Vj {j = 1,. ..,k) are called Kraus operators For any quantum channel acting in an 

TV— dimensional Hilbert space the number of operators k needs not exceed iV^ ^|. Each quantum channel can be 
treated (but not necessarily uniquely) as a pure or mixed states homogenous QIFS. Conversely, for each homogenous 
QIFS, formula ifTTll defines a quantum channel. 

If, additionally, X^j^i VjVj = 1 holds, then A(l/A^) = 1/iV, and the map A is called unital. It is the case if all 
Kraus operators are normal, VjVj = VjVj (j — 1, . . . ,k), however, this condition is not necessary. A unital trace 
preserving CP-map is called bistochastic. An example of a bistochastic channel is given by random external fields | l3^ 
defined by 

k 

p'=A,r(p)^^p, [/.pC/t , (17) 

1=1 

where Ui, i — l,2,...,fc are unitary operators and the vector of non-negative probabilities is normalized, i.e., 
Yli=iPi — 1- The Stinespring-Kraus form ifTfiji can be reproduced setting Vi = ^JplUi. Note that the random 
external field ifTTjl may be regarded as a homogenous QIFS of the first kind (with constant probabilities) with k 
unitary maps Gi{p) = UipUj {i = 1, . . . , k). In particular. Example 7 belongs to this class. In the sequel such QIFSs 
will be called unitary. For a unitary QIFS not only is an invariant state of Ajj , but also the measure (5p, is invariant 
for the operator Pu induced by this QIFS. 

Although a unitary QIFS consists of isometrics, the operator Ajj needs not preserve the standard distances between 
any two mixed states. For the Hilbert-Schmidt metric we have 

Dhs{Au{pi),Au{p2)) <Dhs{pi,P2) ■ (18) 

In fact this statement is true for any bistochastic channels as shown by Uhlmann , but it is false for arbitrary 
CP maps, since the Hilbert-Schmidt metric is not monotone |40|. On the other hand, A(j is a contraction for the 
Bures distance (Riemannian) and the trace distance (not Riemannian), which are monotone and do not grow under 
the action of any CP map Choosing for p2 the maximally mixed state p* = 1/A^, which is invariant with 

respect to Au for any unitary QIFS, we see in particular that the distance of any state pi to p* does not increase 
in time. Similarly, the von Neumann entropy given by H(p) = tr(plnp) for p € Mn does not decrease during the 
time evolution H17|l. On the other hand, the inequality in Eq. Ijl8ll is weak, and in some cases the distance may 
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remain constant. The question, under which conditions this inequality is strong, is related to the problem, for which 
unitary QIFSs the maximally mixed state is a unique invariant state of A^/. This is not the case, if all operators 
Ui commute, since then all density matrices diagonal in the eigenbase of Ui are invariant. Such a situation may 
occur also in subspaces of smaller dimension. To describe such a case we shall call unitary matrices of the same size 
common block-diagonal, if they are block-diagonal in the same basis and with the same blocks. The uniqueness of 
the invariant state of a unitary QIFS is then characterized by the following proposition, the proof of which is provided 
in the appendix A. 

Proposition 2. Let us assume that all probabilities pi (i — I, . . . ,k) are strictly positive. Then the maximally 
mixed state p* is not a unique invariant state for the operator Au if and only if unitary operators Ui (i — 1, . . . , k) 
are common block-diagonal. 

It follows from the proof of this proposition that in this case there exists p ^ p* such that 5p is an invariant measure 
for the operator Pu induced by the QIFS. 

To show an application of Proposition 2 consider a two level quantum system, called qubit, which may be used to 
carry a piece of quantum information. Let us assume it is subjected to a random noise, described by the following 
map: 



Au{p) = (1 - J3)/0 + I [a-ipcri + cr2P<J2 + o^ipa^ 



(19) 



This bistochastic map, defined by the unitary Pauli matrices Ci, is called depolarizing quantum channel [421, and 
the parameter p plays the role of the probability of error. This map transforms any vector inside the Bloch ball 
toward the center, so the length of the polarization vector decreases. In formalism of QIFSs this quantum channel is 
equivalent to the following example. 

Example 9. = V2, k = A, Ui — 1, U2 — cri, U3 — 02, U4, — (T3, pi — 1 — p and P2 — Ps — Pi = P > 0- Since the 
Pauli matrices are not common block-diagonal, the maximally mixed state is a unique invariant state of the CP 
map H19|l associated with this unitary QIFS. 

To introduce an example of QIFS arising from atomic physics, consider a two level atom in a constant magnetic 
field Bz subjected to a sequence of resonant pulses of electromagnetic wave. The length of each wave pulse is equal 
to its period T and it interacts with the atom by the periodic Hamiltonian V{t) — V{t -\- T). Let us assume that 
each pulse occurs randomly with probability p. Thus, the evolution operator transforms any initial pure state by the 
operator 

Ui = cxp{-iHoT/h) (20) 

in the absence of the pulse, or by the operator 



U2 ~ C exp 




(21) 



in the presence of the pulse. The unperturbed Hamiltonian Ho is proportional to B^Jz {Jz is z component of the 
angular momentum operator) and C denotes the chronological operator. Thus, this random system may be described 
by the following QIFS. 

Example 10. fl — V2, k — 2, pi = 1 — p and p2 — p, the Floquet operators Ui Il20|l and U2 (I21|l as specified 
above. The maximally mixed state p* — 1/2, corresponding to the center of the Bloch ball, is the invariant state of 
the Markov operator given by Eq. ifTzjl . For the case of a generic perturbation V, the matrices Ui and U2 are not 
common block-diagonal, and so p* is the unique invariant state for operator II17II related to the QIFS. 

The QIFSs arise in a natural way if considering a quantum system acting on Hn coupled with an ancilla: a state 
in an auxiliary m-dimensional Hilbert space Hm, which describes the environment. Initially, the composite state 
describing the system and the environment is in the product form, tr = p^ ® pf , where pf = Im/m is the maximally 
mixed state, but the global unitary evolution couples two subsystems together. A unitary matrix U of size Nm acting 
on the tensor space Hn ^ Hm may be represented in its Schmidt decomposition form as C/ = J2i^=i y/Qi^i^ ^ 
where the number of terms is determined by the size of the smaller space, K — min{A^^, m^}; the operators and 
act on Hn and Hh respectively, and the Schmidt coefficients are normalized as Yl!i=\ 9» — 1- Restricting our 
attention to the system A one needs to trace out the variables of the environment B which leads to the following 
quantum channel (and to the respective homogenous QIFS): 

K 

p'a - K{pa) = tMUaU^) = 1^^'^PAy^^^ ■ (22) 
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Since for = 1^/^ we have A(p^) = trB(U(p^ ® pf)'U^) = pt, the CP-map A is bistochastic. 

IV. QUANTUM CLASSICAL CORRESPONDENCE 

To investigate various aspects of the semiclassical Umit of the quantum theory it is interesting to compare a given 
discrete classical dynamical system generated by / : 1^ — > 51 with a family of the corresponding quantum maps, 
usually defined as F/v : TCn 'Hn with an integer N . Several alternative methods of quantization of classical maps 
in compact phase space have been applied to construct quantum rnaps corresponding to baker map on the torus 
[4^l43|. Arnold cat map and other automorphisms on the torus periodically kicked top j47| and baker map 
on the sphere [4S] . 

To specify in which manner the classical and the quantum maps are related, it is convenient to introduce a set 
of coherent states \y) £ Hn, indexed by classical points y of the phase space fl. (For more properties of coherent 
states and a general definition consult the book of Perelomov [42|.) They satisfy the resolution of identity formula: 
J^\y){y\dy — 1, and allow us to represent any state p by its Husimi representation, H{y) = {y\p\y)i {y € 51). 
Quantization of a classical map /, which leads to a family of quantum maps Fn is called regular, if for almost 
all classical points x the classical and the quantum images are connected in the sense that the normalized Husimi 
distribution of the state Fj\;\y) integrated over a finite vicinity of the point f{y) tends to unity in the limit N ^ oo 
Another method of finking a classical map with a family of quantum maps is based on the Egorov property, 
which relates the classical and the quantum expectation values [5lll5^. 

In a similar way we may construct QIFSs related to certain classical IFSs. More precisely, a sequence of pure states 
QIFS = {Vn] Fi^N,Pi,N ■ i = l,...,k} induced by two sets of linear maps Vi^N, Wi^N ■ Hn Hn {i = 1, ■ ■ ■ ,k) 
(see © and ©) is a quantization of a classical IFS JFci = {51; Fi,pi : 1, . . . , k}, when: 

• the functions F^jv are quantum maps obtained by quantization of the classical maps fi; 

• the probabilities Pi^N computed at coherent states \y) fulfifi 

P^,Ni\y){y\) = \m^N{\y))f-^P^{y) for y e 51 and z = 1, A: . (23) 

To ifiustrate the procedure let consider random rotations on the sphere, performed along x or z axis. This special 
case of Example 4 may be easily quantized with the help of the components Ji {i = x, y, z) of the angular momentum 
operator J, satisfying the standard commutation relations, [Ji, Jj] — tijkJk- The size of the Hilbert space is determined 
by the quantum number j as iV = 2j + 1 . 

Example 11. k = 2, random rotations are given as the following 

a) classical, J-qi = {51 = S^, fi — Rz(9i), f2 = Rx{02), Pi = P2 = 1/2}- The Lebesgue measure on the sphere is an 
invariant measure of this IFS. 

b) quantum, Tn = {51 = Vn, Fi = exp(i6'i J^), F2 — eyip{i62Jx), Pi = P2 = 1/2}. Since both unitary operators 
are not common-block diagonal, due to Proposition 2, the maximally mixed state p* is a unique invariant state for 
operator lfTT|l to the QIFS J^n- 

A quantization of an IFS of the second kind is given by the following modification of the previous example. 

Example 12. k = 2, random rotations on the sphere with varying probabilities depending on the latitude 6 
computed with respect to the z axis. 

The spaces and the functions are as in Example 11, but 

a) classical IFS J^ci- Pi — + cos0)/2 and p2 = (1 — cos0)/2; 

b) quantum IFS J^n- Pi = 1/2 -I- {Jz)/2j and p2 — 1/2 — {Jz)/2j with N — 2j + 1. Interestingly, this modification 
infiuences the number of invariant states of the IFS. Since p2 vanish at the north pole, = 0, of the classical sphere 
S'^, this point is invariant with respect to Tci- Similarly, the corresponding quantum state \j,j) localized at the pole 
is invariant with respect to the QIFS J^n- 

The above examples of unitary QIFS dealt with simple regular maps — rotations on the sphere. However, an IFS 
may also be constructed out of nonlinear maps, which may lead to deterministic chaotic dynamics. For instance, one 
may consider the map describing periodically kicked top. It consists of a finear rotation with respect to x axis by angle 
a and a nonlinear rotation with respect to z axis by an angle depending on the z component. In a compact notation the 
classical top reads, 101(0;, /?) := Rz{zf3)Rx{a), while its quantum counterpart, acting in the N = 2j + 1— dimensional 
Hilbert space can be defined by TQ{a,[3) :— exp(— i/?J^/2j) exp(— ■jaJj:) [43. This quantum map becomes one of the 
important toy model often studied in research on quantum chaos [^. A certain modification of this model, in which 
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the kicking strength parameter (3 was chosen randomly out of two values, was proposed and investigated by Scharf 
and Sundaram This random system may be put into the QIFSs formaHsm. 

Example 13. Randomly kicked top. 

a) classical, Tc\ = {Vl = S'^, fi = Tci{a,f3), /2 = Tci{a,f3 + A), pi = P2 = 1/2}. 

b) quantum, Tn = {ft = Vn, Fi = TQ(a, /3), /2 = 7q(q!, /3 + A), pi = P2 = 1/2}- For a > and a positive A both 
unitary operators are not block-diagonal, so the maximally mixed state p* is a unique invariant state for operator 
ifTTll related to the unitary QIFS. Our numerical results obtained forQ! = 7r/4,/3 = 2 and A — 0.05 suggest that the 
trajectory of any pure coherent state converges to the equilibrium exponentially fast. 

To discuss a quantum analogue of an IFS with a fractal invariant measure consider the classical IFS presented in 
Example 3. The classical phase space 51 is equivalent to the torus. For pedagogical purpose, let us rename both vari- 
ables x,y into q,p, represented canonically coupled position and momentum. We shall work in TV = 3L— dimensional 
Hilbert space. Let \j)q with j = 1, N be eigenstates of the position operator, and similarly \l)p with I = 1, N be 
the eigenstates of the momentum operator. Both bases are related by \l)p — ^2^=1 where the matrix W is the 

N point discrete Fourier transformation with Wij = (l/-\/]V)e~^^*'-'/^. The classical map /i in Eq. JSJ, representing 
a three-fold contraction in the x direction, corresponds to the transformation Gi of the density operator given by 



Gi{p) = E 1^)9 ( E + m\,p\3j + 72) A {jU ■ (24) 

— l \m,n— / 



In a similar way, the quantum map G2 corresponding to /2 is defined by 

G2{p) = E 1*)^ ( E (3* + ^\iP\^J + ) 0% ■ (25) 

i,j=2L + l \m,n=a / 

The maps G3 and G4 are obtained in analogous way like Gi and G2 , using the eigenstates of momentum operator 

\k)p, 

Gsip) = E 1^=)^ ( E + + n)p\ {l\p , (26) 

k,l=l \m,n=0 ) 

3L / 2 \ 

G4(p)= E |fc)p E (3fc + ™U|3Z + n)J (Zip . (27) 

fe,i=2L+l \m,n=0 / 

The random system defined below may be considered as a QIFS related to the IFS introduced in Example 3. 

Example 14. Quantum tartan specified by the following QIFS: Tt^ — — Vn , k — A, Gi, G2, G3, G4; pi — P2 = 
P3=Pi^ 1/4}. 

An invariant states for the maps A induced by this QIFS are illustrated in Fig. 1 for iV = 3"^, = 3^ and A^ = 3^. 
Invariant quantum state p, is shown in the generalized Husimi representation 



HM-^^-f^, (28) 
27r {q,p\q,p) 

based on the set of coherent states on the torus \q,p) = Y^p^^/^X^'i~^/^\k) . The reference state \k) is chosen as 
an arbitrary state localized in (1/2, 1/2) 

{n\K) = (2/Ar)-i/4g-^(«-A'/2)^/A'-z7™ ^ ^29) 

while X denotes the operators of shift in position X\ j) = \j + 1) , with an identification \ j + N) ~ \ j) for j = 1, . . . , A^. 
Similarly Y shifts the momentum eigenstates, Y\l) = \l + 1) and \l + N) = \l) for I = 1, . . . ,N. The quantum state 
\q,p) is well locaHzed in the vicinity of the classical point (g,p) on the torus This representation of quantum 

states corresponding to the classical system on the torus was used in the analysis of an irreversible quantum baker 
map | i56j . 
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Figure 1: "Tartan-like" invariant density of the QIFS defined in Example 14 for (a)A'' = 3"*, (b) N = 3^, and (c) N = 3^- 
dimensional Hilbert space, shown in the generalized Husimi representation. Invariant measure of the corresponding classical 
IPS on the torus Eq. I^J occupies a fractal set (d). 

The larger value of N, the finer structure of the invariant state is visible in the phase space. In the semiclassical 
limit N ^ oo, (which means h ^ 0) the invariant state p* tends to be locaHzed at the fractal support of the 
invariant measure of the classical IFS, shown for comparison in Fig. Ic. Strictly speaking, for any finite N, the 
Husimi distribution of the quantum state p, does not posses fractal character, since self-similarity has to terminate 
at the length scale comparable with ^/h. In other words, quantum effects are responsible for smearing out the fractal 
structure of the classical invariant measure. However, the classical fractal structures may be approximated with an 
arbitrary accuracy by quantum objects in the semiclassical limit [H^ . 

V. CLOSING REMARKS 

Classical iterated function systems display several interesting mathematical properties and may be applied in various 
problems from different branches of physics. In this work we have generalized the formalism of IFSs introducing the 
concept of QIFSs. Quantum iterated function systems may be defined in the space of pure states on a finite dimensional 
Hilbert space Hn, or more generally, in the space of density operators acting on Hn- As their classical analogues, 
QIFSs allow a certain degree of stochasticity, in the sense that at each step of time evolution the choice of one of the 
prescribed quantum maps is random. 

This formalism is useful to describe several problems of quantum mechanics, including non-unitary dynamics, 
processes of decoherence and quantum measurements. In fact, the large class of quantum channels, called random 
external fields may serve directly as examples of a QIFS. Furthermore, for several classical IFSs one may construct the 
corresponding QIFSs and analyze the similarities and differences between them. As shown in the last example, one 
may focus on the fractal properties of invariant measures of some classical IFSs and study their quantum counterpart. 
Thus the concept of QIFS allows one to investigate the semiclassical limit of random quantum systems. 
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Appendix A: PROOF OF PROPOSITION 3 

We start from the following lemma: 

Lemma 1. Let U = {Unm)nm=i n be an A/'— dimensional unitary matrix. Assume that there exist two non- 
empty sets of indices A and B such that: A(J B = I := {1, . . . ,N} and AnB = Then, Unm = for n e A and 
m £ B, implies Unm = for n e _B and m ^ A. 

Proof of the lemma: We compute the number of elements of the set A: 

■'nm I 



1^1 = E E 
= EEit^"'"i' + EEit^' 

uGA mGA 

= EEif^' 



neA mel 

2 

nm I 

nGAmGA neAmeB 
2 



nm I 

nGA mGA 

EE if^"-i'-EEi^"-i' 

nel meA neB meA 

1^1 - E E 

neB meA 



,2 



and so J2neB J2meA Pnm.l'^ = 0, as required. 

Now we turn to the proof of Proposition 3. 

Let Ui (?' = 1, , .... A;) be block-diagonal in the common base, and let dimension of the blocks be ai, . . . ,aL, 
where J2^=i — ^- Define a diagonal density matrix as a direct sum 



■ 1 '^3 



where X^^^i ctj — 1- Then, UipUj = p for every i = l,...,k. Hence p is A[/— invariant and 6p is a P;/— invariant 
measure on Vn for an arbitrary choice of {crj)j—i j^- 

<=) Let p be an invariant state for h.u such that p^ p^. Then p can be written in the form 

JV 

p = ^a„|*„)(*„| , (A2) 



where e Vn, (*n|*m) = 5nm {n,m = 1,...,N), and cti < C72 < ••• < ctjv; (Ti < 1/N. For 7 e [0,1] the 

density operator p' = 7/3 + (1 — 7)^, = J2n=i ^nl^n)(^n|) where tr^ = 7(t„ + (1 — 7) (n = 1, .... TV) is also an 
invariant state for Ay. Put 7 := 1/(1 — ctiA'^). This choice implies a[ = and J2n=i ^ri = 1- Assume that cr^ = for 
n = 1, . . . , n' and a'^ > for n = n' + 1, . . . , N , where n' > 1. The equation Au{p') = p' can be rewritten in the form 



k N 



'^n = EPi E K^'WI"'^™ ' (A3) 



i—1 m—1 



where {Ui)nm {n,m = 1, . . . , N) are the elements of matrices Ui {i = 1, ... ,k) in the basis (|^'n))„=i jv- 
For n = 1, . . . , n' we get 



N 



^ = J2p^ E mi)nm\''cT'^ . (A4) 
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Hence {Ui)nm = for n = l,...,n' and m = n' + I, . . . , N . Using Lemma 1, we deduce that {Ui)nm = for 
n = n' + 1, . . . , iV and to = 1, . . . , n'. Thus Ui {i — 1, . . . ,k) are common block-diagonal. 
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